Fluxes of cosmic rays: A delicately balanced stationary state 
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; Abstract 

m 

£NJ . The analysis of cosmic rays fluxes as a function of energy reveals a knee 

. slightly below 10 16 eV and an ankle close to 10 19 eV. Their physical origins 

. remain up to now quite enigmatic; in particular, no elementary process is 

o . 

CN ■ known which occurs at energies close to 10 16 eV. We propose a phenomeno- 

O 

logical approach along the lines of nonextensive statistical mechanics, a for- 
q . malism which contains Boltzmann-Gibbs statistical mechanics as a particular 

h'. 

case. The knee then appears as a crossover between two fractal-like ther- 

o3 



X 



mal regimes, the crossover being caused by process occurring at energies ten 
million times lower than that of the knee, in the region of the quark hadron 
transition (~ 10 9 eV). This opens the door to an unexpected standpoint for 
further clarifying the phenomenon. 

PACS number: 95.30.-k, 95.35.+d, 95.85.Ry, 05.90.+m 

Cosmic rays fascinate since long. They provide galactic, extragalactic and cosmological 
information, related to recent or very old events concerning various sources, going back to the 
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early times of the universe [1]. They reflect all types of elementary process and interactions, 
and are associated with phenomena of very different space and time scales. The complete 
physical scenario still remains quite enigmatic, although a variety of specific mechanisms for 
acceleration and propagation have been advanced along the years for various energy regions. 
The most known of these mechanisms is the Fermi one [2], which addresses acceleration in 
magnetized turbulent plasma, but many others have been advanced [3-10] in relation with 
the knee and energies below it; for the energies beyond these, for the ankle, as well as for 
general reviews, see [11-13]. 

Through various types of detectors, the flux of cosmic rays at the top of the Earth at- 
mosphere has been measured [14] and varies from 10 4 down to 10~ 29 [m 2 sr s GeV] -1 for 
energies increasing from 10 8 up to near 10 21 eV: See Fig. 1. This distribution (which spans 
13 decades in energy and 33 decades in flux!) is not exponential, hence it does not corre- 
spond to Boltzmann-Gibbs (BG) statistical mechanics thermal equilibrium. Consistently, 
even at a phenomenological level, i.e., without specifying any concrete model or mechanism, 
this problem represents a challenge. This is the one we address here. We shall use a point of 
view based on a current generalization of Boltzmann-Gibbs statistical mechanics, referred 
to as nonextensive statistical mechanics, we shall briefly describe later on. The first step 
will be to remark that the fluxes of cosmic rays in general, and the studies of the "knee" and 
the "ankle" in particular, involve phenomena such as turbulence (see, for instance, [15]), 
anomalous diffusion and fractality (see, for instance, [16]), self-organized criticality (see, for 
instance, [17]), long-range interactions (classical and quantum gravitation), among other 
complex phenomena (such as, for example, possible nonmarkovianity [18]). It is precisely 
such phenomena that constitute the scope of nonextensive statistical mechanics; for turbu- 
lence and related matters see, for instance, [19-21]; for anomalous diffusion and fractality 
see, for instance, [22-28]; for self-organized criticality see, for instance, [29]; for long-range 
interactions see, for instance, [30-32]. 

Let us now first briefly review the usual, BG thermostatistics. If we optimize under 
appropriate constraints the BG S = —k^^ilnpi {k = Boltzmann constant; {pi} = micro- 
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scopic probabilities) we obtain the celebrated equilibrium distribution p t = 6 z oc e~~^ E% 
(P = 1/kT, Ei = energy of the i-th state; Z = J^ . = partition function). Excepting 

for the trivial normalizing factor 1/Z, this distribution can alternatively be obtained as the 
solution of the linear differential equation 

dpi/dEi = -/3pi . (1) 

In order to deal with a variety of thermodynamically anomalous systems, a more general 
formalism, nonextensive statistical mechanics, was introduced in 1988 [33-36]. It is based on 
the generalized entropic form S q = k(l — J2iPi)/(q — 1) (Q e ^ an d <Si = S). Its optimization 
under appropriate constraints yields [37] a power-law, pi oc [1 — (1 — q)/3 q E i \ 1 -i = e q q 1 
(definition), which recovers the BG weight for q — 1 (f3i = 0). As usual, kT q = 
characterizes the conveniently averaged energy. This anomalous equilibrium-like distribution 
can be alternatively obtained (excepting for the normalizing factor) by solving the nonlinear 
differential equation 

dpi/dEi = -(3 q p1 . (2) 

This generalized weight naturally emerges in ubiquitous problems such as fully developed 
turbulence [19-21] (which is relevant for the Fermi mechanism, and most probably for others 
as well), electron-positron annihilation [18], motion of Hydra viridissima [38], long-range 
many-body Hamiltonians [31], among many others. 

We now use the above differential equation path in order to further generalize the anoma- 
lous equilibrium distribution, in such a way as to have a crossover from anomalous (q ^ 1) 
to normal (q = 1) thermostatistics, while increasing the energy. We consider then the 
differential equation 

dpi/dEi = -p lPi - (/3 q - 0)p\ , (3) 

whose solution is pi oc [1 — |^ + |^ eO?- 1 )/ 3 !^-^!, The crossover typically occurs for q > 1 
and Pi P q , the distribution being anomalous at low energies and BG at high energies. 
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It is undoubtedly interesting to notice that this differential equation precisely coincides, for 
q — 2, with the heuristic one that in 1900 led Planck to the discovery of the black-body 
radiation law and ultimately to quantum mechanics [39]. 

Finally, by doing one more step along the same direction, we can further generalize the 
differential equation, now becoming 

dpi/dEi = -Ptfpf - (fig - . (4) 

This manner of writing the coefficient of the p\ term has the advantage of recovering the 
simplest generalization of the BG distribution by considering fi q i = P q or p q > = or even 
q' = q. For 1 < q' < q and P q i <C P q , a crossover occurs at 

E crossover = [(q - l)P q 0/[(q' - l)P ql \^ . (5) 

For E <C E crossover we have an anomalous distribution characterized by (q, P q ) (namely 
Pi oc e q q l ), whereas for E ^> E crossover we have a different anomalous distribution char- 
acterized by (q',P q >) (namely pi oc e ql q ' l ). The exact solution of the above differential 
equation (the most general one considered here) is given by Pi oc f(Ej) where / _1 (x) is an 
explicit monotonic function of x involving hypergeometric functions (see Ref. [40] for details). 
Interestingly enough, this precise solution arrives in the discussion of the re-association of 
CO molecules in Myoglobin [40], where time plays a role very analogous to the one played 
by energy in our cosmic rays problem. This time-energy analogy is not surprising after all 
if we take into account that, in the history of the universe after the big-bang, the time scale 
reflects the energy scale, as discussed in detail in Ref. [41]. 

The flux $(£■) can be obtained straightforwardly from pi oc f(Ej) by calculating the 
density of states uo(E). In the ultrarelativistic limit E oc |p| (p being the momentum), 
which we adopt here for simplicity given the high values of the involved energies, the density 
of states of an ideal gas in three dimensions is given by uo(E) oc E 2 , hence $(E) = AE 2 f(E), 
where A is a normalizing factor (and where red shift effects have been neglected). With this 
expression we fit the observational data and obtain the results displayed in Fig. 1. As we 
can see, the agreement is quite remarkable. 
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Our summarizing comments are: 

(i) The high quality agreement over so many decades, including crossovers between dif- 
ferent regimes, suggests that the phenomeno logical approach is correct, and specific models 
clarifying the various physical mechanisms that are involved should essentially satisfy it; 

(ii) The deep explanation of the knee might well be found at energies extremely lower 
(ten million times lower, in fact), basically at energies related to the characteristic tempera- 
tures obtained from the fitting, namely 9.615 10 7 eV (energy comparable to the pion mass) 
and 1.562 10 9 eV (energy corresponding to the quark-hadron transition [41], as well as to 
the proton mass), and given the fact that q and q' differ by only 3% (which in log-log rep- 
resentation asymptotically determines two almost parallel straight lines which very slowly 
approach to each other and intersect at E = E crossover ). The existence of two thermostatis- 
tical regimes, respectively related to (q,/3 q ) and (q',j3 q >), could correspond to two different 
mechanisms of acceleration/propagation, for instance related to galactic and extra-galactic 
contributions. It is worthy stressing at this point that the location of the knee emerges here 
through Eq. (5), i.e., using only the basic four phenomenological parameters (q, (3 q , q', (3 q i), 
and not by introducing an extra parameter on top of the previous ones (whose role would 
be to fix the knee at the right value); 

(iii) Since the entropic index q is known to reflect (multi) fractality [42,43], the present 
results strongly suggest that either the generation or the transport (or both) of cosmic rays 
occur in scale invariant media, which is consistent with Ref. [16]; 

(iv) With the better observational statistics (i.e., with refined precision) at these very 
high energies, expected from the Pierre Auger Observatory (or from similar projects), it 
might happen that the ankle disappears; if it does not, then it is probably associated with 
strongly nonstationary phenomena (perhaps related to memory effects from early cosmic 
stages), certainly out from the present thermodynamical description — such study would 
eventually clarify the physical meaning, if any, of the GZK feature [44-47] (see also [49]). 

(v) The Kaskade collaboration [9] has shown the relevance of the various constituents 
which compose cosmic rays. In order to take this into account, the present approach could 
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be further improved by using the grand-canonical instead of the canonical ensemble used 
here. 

(vi) The average energy (E) = J °° dE E $(£')/ / °° dE $(£') has been calculated, within 
the present approach, to be (i?) ~ 2.489 GeV. Any connection of this value with other 
cosmological or astrophysical quantities is of course very welcome. 

Stimulating remarks from I. Bediaga, L. Galgani, R.C. Shellard and L. Masperi are 
acknowledged. Partial support from CNPq, PRONEX, CAPES and FAPERJ (Brazilian 
agencies) is also acknowledged. 
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FIGURES 
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FIG. 1. Energy dependence of the fluxes of cosmic rays. Experimental error bars are indicated 
whenever available. The continuous curve is the one we obtain within the present phenomenological 
approach. The dashed curve is an optimized BG one (even at relatively low energies it fails by 
very many decades). The knee corresponds to E crossover . Inset: Linear-linear representation of the 
low energy fluxes. 
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